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(NJ . ABSTRACT 

' We study the universality in gravitational waves emitted from non-rotating neutron 

D , stars characterized by different equations of state (EOS). We find that the quasi- 

£^ ' normal mode frequencies of such waves, including the w-modes and the /-mode, dis- 

play similar universal scaling behaviours that hold for most EOS. Such behaviours are 
shown to stem from the mathematical structure of the axial and the polar gravitational 
wave equations, and the fact that the mass distribution function can be approximated 
^ ' by a cubic-quintic polynomial in radius. As a benchmark for other realistic neutron 

stars, a simple model of neutron stars is adopted here to reproduce the pulsation 
frequencies and the generic scaling behaviours mentioned above with good accuracy. 
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1 INTRODUCTION 

The search for gravitational wave has been the goal of endeavour for generations of physic ists since its existence was predicted 
in the theory of general relativity (see, e.g. a review by Thorne in lHawking fc Israeli Il987i chap. 9). In the beginning of a new 
millennium, the breakthroughs in technology and endless efforts of researchers strongly convince us tha t the first detection of 
gravi tational waves emitted from violent stellar collapses will be realized within several years (see e.g. iHugheslEool iMasonl 
2004, and references therein). Gravitational wave detectors of different designs, including ground-based interferometers (e.g. 
LIGO and VI RGO), and resonant antenna s (e.g. EXPLOPER and NIBOE), are either already in operation or in their final 
testing stages jHughesJbooal lMasonll2004L These detectors, together with the space-based interferometers LISA, which is 
currently under construction and is expected to be launched around 2011, can detect gravitational waves in different frequency 
ranges. As a consequence, gravitational-wave astronomy that enables researchers to pro be astronomical activities at extremely 
large distances from gravitational waves detected emerges as a new subject of interest faughedkool lMasonll200i) . 

Among all the possible source s of gravitational waves, pulsating neutron stars may probably be the most interesting 
one iThorne fc Campolattarolll967l) . It is believed that gravitational waves emitted in binary neutron mergers are plausibly 
detectable with the current tech nology and the frequency of detectio n could be as high as several hundreds per year after 
upgrading the present detectors jHugheJbnoal IBelczvnski et alJbOQlh . Besides, waves generated in the formation of neutron 
stars via stellar core collapses might als o be strong enough for detection provided that the collapse is sufficiently asymmetric 
iFrver et alJboollLindblom' et al.ll998l) . Most interestingly, such gravitational waves are likely to carry the information about 
the internal structure of neutron stars from which they are emitted. As the structure of a neutron star, such as its radius and 
density distributio n, certainly reflects the properties of nuclear and quark matters, a close exami nation of the wave signals is 
deemed rewarding l|Andersson fc Kokkotaslll99rl Il99sl. iBenhar et alJll999l. IrCokkotas et alJbOQlh . 

Gravitational waves are commonly analyzed in terms of quasi- normal modes (QNMs), which properly describe damped 
harmonic pulsations. A QNM oscillation has a time dependence exp(wt) and is c haracterized by a complex eigenfrequency ui — 
uJr + iuJi l)Presslll97ll . lLeaveilll98(ilChing et alJll99rllKokkotas fc Schmidtll999l) . The QNM frequencies of gravitational waves 
of a neutron star with a fixed mass M axe generally model-sensitive and significantly depend on the EOS adopted in the stellar 
model. However, some universal behaviours in the frequency w r and the damping time r = l/u)j of the leading gravitationa l 
wave ui-mode of non-rotating neutron stars have recently been observed ijAndersson fc Kokkotasl|l998L IBenhar et ai] ll999V 
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which can be summarized by the following pair of formulas: 
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worthy of remark that such simple scaling formulas ap ply to both axial and polar gravitational wave io-modes of neutron 



Here R is the radius of the star, a r , fo r , ai, bi and c\ are model-independent constants determined from curve fitting. It is 

h y of remark that such simple scaling 
stars lAndersson fc Kokkota3ll99gL iBenhar et alJll999 

On t he other hand, QNM frequenc ies of the fluid /-mode are discovered to follow approximately another pair of universal 
formulas lAndersson & Kokkotas 1998): 
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where again a T , ft T , en, and fti are model-independent constants chosen to yield best fit. Based on the above-mentioned 
universal behaviours, it has been shown that the radius and the mass of a neutron star can be inferred from the pulsation 
frequencies of its fundamenta l fluid /-mode and t he first w-mode, and the equation of state (EOS) can in turn be identified 
lAndersson fc Kokkotaslll998l iBenhar et aflll999h . 

These universal behaviours are intriguing in their own right. Why can Eqs. successfully reproduce the QNM 

frequencies for stars with different EOS? What are their physical interpretations? Is there any exception to them? In this 
paper, we seek the physical mechanism underlying the above-mentioned universality in the gravitational wave of non-rotating 
neutron stars. By properly scaling the axial and polar wave equations, we show that the scaled complex eigenfrequencies 
Mui of axial and polar w-modes, and /-mode oscillations, to a good approximation, depend only on the compactness M/R. 
Moreover, to interpret such EOS-independent generic behaviour, we use a si mple stellar m odel — the Tolman VII model 
(TVIIM) — to approximate the mass distribution inside realistic neutron stars jTolmaiJl939h . and find that TVIIM is indeed 
a good approximation. We then look for the QNMs of TVIIM with the scaled wave equation. Interestingly enough, such 
QNMs manifestly depend only on the compactness of the star and, in addition, demonstrate the scaling behaviour discovered 
in realistic neutron stars. Therefore, the universality in the QNMs (™-mode and /-mode) is ascribable to the fact that that 
the mass distribution of most realistic neutron stars can be well approximated by TVIIM. 

The organization of our paper is as follows. In Sect. 2, we show how the frequencies (including the real and the imaginary 
parts) of QNMs of neutron stars with different EOS can be captured by a single formula. In Sect. 3 we briefly review the 
equilibrium configuration of a neutron star and introduce TVIIM as an approximation of realistic neutron stars. Sections 4 
and 5 study the scaling behaviours of the axial and the polar pulsation modes, respectively. In Sect. 6 we discuss potential 
application of our finding and consider some possible exceptional cases where obvious deviations from the universality are 
expected. We then conclude our paper in Sect. 7. Unless otherwise stated, geometrized units in which G = c = 1 are adopted 
in the following discussion. 



2 GENERIC BEHAVIOUR OF QNMS 

The th eory of non-radial oscillations of relati vistic stars was pioneered bv lThorne fc Campolattarol l)l967h . Later. lLindblom fc Detweilerl 
ll983l) . and lChandrasekhar fc Ferraril Jl99lh further simplified the equations describing gravitational wave generation and fluid 
motion, which will be detailed in the later part of this paper. To display the universality in the oscillation modes, we numeri- 
cally integrate the gravitational wave equations for neutron stars, and impose the outgoing wave boun dary con dition at sp atial 
infinity to locate relevant QNMs. Eight different EOS, including models A lPandha ripandelll971 al) and C l|Pandharipandel 
Il971bh proposed b y Pandharipande, t hree models (AU, U U and UT) propos ed by I Wiringaetal.1 (|l988l) . models APRl and 
APR2 proposed bv Ukmal et al.1 l|l998h . and model GM24 l)oiendenninJl997l . p. 244) are considered in our papers. Figures [Tl 
[5] and [3] show the relation between Mu and the compactness M/R for the first axial ui-mode, the first polar w-mode and the 
/-mode, respectively. It is clearly shown that both Mui r and MuJi depend on the compactness in a universal way. In particular, 
as shown by the solid line in these figures, a quadratic function in M/R, namely: 

'M\ 2 , /M s 



— (f) -'(f)- 

successfully reproduces the dependence on the compactness in each case. Here a, b and c are complex constants determined by 
regression and their values for individual cases are tabulated in Table 1 for reference. While reproducin g the observed universal- 
ity of the wi-mode as described in Q and Eq. JHJ is different from © and Q discovered previously (lAndersson fc Kokkotasl 
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Table 1. Values of a, b and c in l5l . 



Wave type 


Data 




a 


b 


c 


1st axial m-mode 


realistic stars 


-3.9 


-i5.6 


2.8+il.6 


-0.030+i0.12 




TVIIM 


-4.4 


-16.3 


3.1+il.9 


-0.072+i0.098 


2nd axial to-mode 


realistic stars 


-13 


-iS.l 


6.1+i2.8 


0.034+i0.11 




TVIIM 


-14 


-i9.3 


6.7+i3.3 


-0.049+i0.060 


1st polar ui-mode 


realistic stars 


-8.2 


-ilO 


3.9+i3.3 


0.055+i0.043 




TVIIM 


-8.6 


-ill 


4.2+i3.6 


0.0080+i0.014 


2nd polar ui-mode 


realistic stars 


-18 


-i8.6 


8.0+i3.1 


0.031+i0.11 




TVIIM 


-20 


-i9.9 


8.8+i3.6 


-0.077+i0.058 


polar /-mode 


realistic stars 


0.15 


-i5.8E-4 


0.56+i6.7E-4 


-0.020-i6.2E-5 




TVIIM 


0.14 


-16.9 E-4 


0.60+i7.2E-4 


-0.027-i6.8E-5 



1998). In fact, for nonzero /3 r , Eq. leads to explicit mass-dependence of Mu t , which is not observed in our study. On the 
other hand, we note that p-mode QNM frequencies (not shown in the present paper), especially the imaginary part, do not 
demonstrate analogous behaviour. However, as shown in Fig. 0] and such universality also exists in non-leading axial and 
polar ui-modes. 

The neutron stars considered here are constructed from nuclear matters with different compositions (e.g. pure neutron 
matter and hyperon matter) and different stiffness. It is well known that the mass-radius curve for neutron stars reveals 
marked EOS dependence. Yet the scaled QNM eigenfrequency, Muj, can display an EOS-independent generic behaviour. Two 
questions arise naturally from our finding. First, why is the scaled eigenfrequency, Mui, EOS-insensitive? Second, does the 
universal curve in Eq. J5J has any physical meaning? In the ensuing discussion we will work out the answers to these questions. 



3 EQUILIBRIUM CONFIGURATION AND TVIIM 

To look for a simple yet robust model for neutron stars, we first briefly review the equilibrium configuration of non-rotating 
neutron stars. In spherical coordinates, (t,r,8,ip), the geometry of spacetime around a non-rotating neutron star is given by 
the line element: 

ds 2 = -e" (r) dt 2 + e A(r) dr 2 + r 2 (d9 2 + sin 2 6dtp 2 ). (6) 

where the metric coefficient e A ' r ' is determined by the mass distribution function m(r), the mass inside circumferential radius 
r, as follows: 

e . Mr) = : _2m(r)^ 
r 

On the other hand, the metric coefficie nt e 1 ^ 1 "^, the density of mass-energy, p(r), a nd the pressure P(r) obey the Tolman- 
Oppenheimer-Volkoff (TOV) equations l|Tolmanlll939l . lOppenheimer fc Volkofdll93flh : 

du 2m + 87rr 3 P 

dr r(r — 2m) ' 

dm 2 /„n 

— = Airrp, (9) 

f - ('») 

These four equations, together with an EOS, P = P{p), and the following boundary conditions at the surface of the star 
where r = R: 

P(r = R)=0, (11) 

e"*W = e^> = 1 - (12) 

suffice to determine the equilibrium configuration of a neutron star of mass M and the spacetime around it as well. 

For a given physical EOS and a given central density p(r = 0) = po, this system of equations can be solved by numerical 
integration. However, the probl em can also tackled by expanding relevant quantities as power series of r about the origin 
llChandrasekhar fc Ferrarilll99ll) : 

p(r) = po + P2r 2 + ... , (13) 

m(r) = 4.(^ + ^1 + ...), (14) 
P(r) = Po + P2r 2 + ... . (15) 
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Similarly, the series expansions of the metric coefficients are readily obtainable from (J7J and JSJl. Direct substitution of these 
expansion i nto the TOV equations and the EOS yields formulas expressing all the expansion coefficients in terms of the central 
density po llChandrasekhar fc Ferrar3ll99ll) . 

Motivated by this expansion scheme that works nicely about the origin r = 0, we consider a mass distribution: 



m c (r) = M 



5/r\ 3 3 / r \ l 
2\r) ~2\Rj 



(16) 



which was first introduced bv lTolmanl (|l939F ) and is commonly referred to as Tolman VII model (TVIIM). In addition to being 
a good approximation near the origin, TVIIM satisfies the boundary condition p(r — R) = at the stellar surface, which is 
correct for most EOS as P(r = R) = 0. Therefore, it is reasonable to conjecture that TVIIM is able to provide a good global 
approximation to the exact mass distribution. The conjecture is indeed confirmed by the results shown in Fig. E] where the 
normalized mass distribution function m(r)/M is plotted against the scaled radius r/R for stars described by different EOS 
with a common compactness M/R = 0.2. It is manifest that m( r) obtained from stars with different EOS are close to that of 
TVIIM, and similar conclusion has also been drawn recently bv lLattimer fc Prakashl J200lh . This discovery suggests that the 
QNMs of TVIIM could be close to those of realistic neutron stars. 

Moreover, TVIIM has a distinct feature that gives rise to the universal behaviours observed in the w-mode QNMs. 
Introducing the scale transformations: 

f = r/M, (17) 

m(f) = m(r)/M, (18) 

P(f) = M 2 P{r), (19) 

p{f) = M 2 p{r), (20) 

*>(r) = Kr), (21) 
we rewrite the TOV equations as follows: 



dm 
df 

di> _ 2m + 8nr 3 P 
df f (f — 2m) 



Amfp, (22) 

(23) 



AP 1, ~^Av 

df = -2^ +P )df- < M > 
These scaled TOV equations have exactly the same forms as the original ones. In particular, for TVIIM the scaled mass 
distribution function is given by: 

m c (r) =m c (r)/M= - J r 3 - - J r 5 . (25) 

It is clear that the function m c (f) is solely characterized by the compactness M/R and is independent of the mass M. 
Consequently, other scaled functions (e.g. /5(f), P(f) and S(f)) associated with TVIIM depend only on M/R as well. In fact, 
this statement holds as long as m(f ) does not contain any explicit mass- dependence. However, as shown in Fig. E] TVIIM 
provi des the simplest, yet robust, global approximation to m(r). Besides, TOV equation for TVIIM is exactly solvable ijTolmanl 
1939), thus further simplifying relevant calculations. In the subsequent discussion, we will show that both the axial and the 
polar gravitational wave equations also exhibit similar scaling behaviour, thus leading to the universality in the QNMs. 



4 SCALING BEHAVIOUR OF AXIAL QNMS 

The equation of motion for relativistic stars can be decomposed in terms of spherical harmonics and in turn classified into 
axial and polar oscillations l|Thorne fc CamDolattarc1ll967i iLindblom fc DetweileJl98.4 Ichandrasekhar fc Ferrarllll99lh. Th e 
equation for the axial oscillations of neutron stars is given by a Regge- Wheeler- type eauatioli llciiariflras^k^^ : 

<|J ' ' (26) 



dr 



T +Lo 2 -V{r t ) 



V(r») =0, 

where the tortoise coordinate r„ and the potential V are defined by 
, = / e ( "^ +A)/2 dr 



(27) 



and 
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V(r») = = 1 [l(l + l)r + 4irr 3 (p-P)-6m(r)}, (28) 

respectively. Outside the star, the pressure and the density vanish and hence the tortoise radial coordinate there can be 
simplified to 

r 



r + 2Mln ("2M 



1 )+C, 



(29) 



where C is a constant that can be obt ained by matching I27t with 1291 1 at r = R. Similarly, the potential V is reduced to the 
well known Regge- Wheeler potential jReeee fc WheeleiHl957f) . 



V rw (r.) 



1(1 + 1) _ 6M 

iy 2 ty* 3 



The axial mode wave equation can be scaled by multiplying M 2 to yielding 



[d* 2 

dff + " 


-V(f t ) &(?„) = o, 


where 




Cu — 


Mlo, 




r*/M, 


■ = 


^(r„), 


and 




W.) = | 


e^f" 3 [l(l + l)f + 4irr 3 (p 
(l-2f- 1 ) [I(i + l)f- a _ 



6f 



- 6m(f )] , 



f sC R/M; 
f > P/M. 



(30) 

(31) 

(32) 
(33) 
(34) 

(35) 



For TVIIM, rh(r) = m c (r), p and P are universal functions of f with the compactness being the unique parameter. As a 
consequence, QNMs of TVIIM depend only on the compactness. For other realistic stars whose mass profiles are similar to that 
of TVIIM, it is plausible to expect that the scaled QNM frequency ci) is close to that of TVIIM with the same compactness. 
Hence, we argue that ui of realistic stars has a trend similar to that of TVIIM. This conjecture is confirmed in Fig. Q where 
QNMs of TVIIM (solid-circles) are shown to be close to those of other realistic stars. It is amazing that the dotted line, the 
best quadratic fit to QNMs of TVIIM, almost coincides with the solid line that is obtained by fitting a quadratic expression 
in M/R to Q for stars with eight different EOS. Obviously, the universal behaviour in QNMs of neutron stars, including 
realistic ones and TVIIM, is well approximated by the dotted line. For reference and comparison, we record the values of the 
regression parameter a, b and c (see of the dotted line in Table 1. One can see that the parameters a and b of the two 
lines are indeed quite close, which agrees with our observation from Fig. On the other hand, there is a disparity between 
the parameters c for the two cases. However, such a disparity is obviously unimportant and does not affect the agreement 
between the two curves. 



5 SCALING BEHAVIOUR OF POLAR QNMS 



The equation of motion for polar case allows pulsating modes of various kinds, including the fundame ntal fluid /-mode, the 
gravi tational wave w-mode, the pressure p-mode and the gravity g-mode (see, e.g, a recent review bv lKokkotas Sz Schmidt! 
1999). To investigate why QNMs of f- mode a nd po lar w-mode oscillations display generic behaviour as discussed above, we 
adopt the notation of iLindblom fc Detweileii l)l983h and consider the set of polar wave equations governing the variation of 
the thre e metric perturbation functio ns, Ho, -Hi, and K, and the two fluid perturbation functions, W and V, inside the star 
(r < R) jLindblom fc Detweileilll983l) : 



K' 
W 
X' 



-Ho + 1-1(1 + l)#i 

r 2r 



(i + X) 



Hi + -e x [H + K~ 16ty( p + P)V], 
r 

8n(p + P)e x/2 



K 



-W, 



(Z + X) 



W + re 



A/2 



-i//2 



7 P 



-X 



1(1 + 1) 



X + (p + P)e 



v/2 



H + 



1 



Ho + K 



At i r>\ A/2 i 2 A/2— v 

4tt(p + P)e ' + oj e ' 



1 2 
2 r 



2 _„ 1 1(1 + 1) 

rtu c + 

2 r 



Hi + 



1 /3 , 

2 



K 



1 1(1- 
2 j 



-v V 



(36) 
(37) 
(38) 

(39) 
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Here 7 is the adiabatic index: 
and X is defined as: 



7 



(40) 



X = uj 2 {p + P)e- v/2 V - -P'e (l " X)/2 W + l(p + P)( 



1//2 



Ho, 



2 ' ^ 

which is related to the Lagrangian perturbation in pressure, Ap, through the equation X — —e v ^ 2 Ap jLindblom fc Detweileil 
Il983h . These four first order ODE's, together with the constraint equation, completely determine the dynamics inside the 
star. 

Outside the star, the whole problem could be described by only two metric perturbation functions, namely Ho an d K. A 
linea r combination of these two functions defines the Zerilli wave function Z(r t ) that satisfies the Zerilli wave equation JZerillil 
1970): 



dr 2 



1 * 
+ CO 



Vz(r*) 



Z(r.) = 0, 



with Vz(r«) being the Zerilli potential 
1 - 2M/r 



V z (r») = 



[2n 2 (n + l)r 3 + 6n 2 Mr 2 + 19,nM 2 r + 18M 3 ] 



(42) 



(43) 



r 3 (nr + 3M) 2 
and n={l- 1)(/ + 2)/2. 

A close examination of the polar wave equations inside the star reveals that the scaling behaviour is likely to be marred. 
The culprit is the presence of the adiabatic index 7, which directly depends on P (instead of P) and is hence not amenable 
to scaling. In fact, it has been shown QNMs of p-mode oscillations (especially the damping rate) do not follow the universal 
behaviour. However, as X is proportional to the Lagrangian perturbation in the pressure, it should be negligibly small for 
ui-mode and /-mode waves. To verify this issue, we solved the polar wave equations with the assumption X = and found 
that eigenfrequencies of ui-mode and /-mode waves are almost unaffected by such assumption (see Fig. for the result in the 
ui-mode case). This shows that the term X/7 in 1381 can be omitted as far as tu-mode and /-mode waves are concerned. 

Under the approximation X = 0, which is valid for w-mode and /-mode oscillations, we can show that the universality 
observed in the axial case also exists in the polar case. Again we have to introduce more scale transformations: 



W/M 2 



W = 

X = M'X, 
V = V/M 2 , 



(44) 
(45) 
(46) 



while for Ho, Hi, K, X, the transformed quantities are identical to the original ones. In terms of these scaled variables, 
Eqs. lOZJ and (jjg) then become 



H[ = -i 

r 



/ + i + ^_ +47r? V(p-p) 



K' = Ih + ±-J(l + l)Hi - 
r 2r 



(l + l) 



Hi + — e [Ho -\- K — 16tt(p + P)V], 
r 

R _ 8n(p + P)e^ 2 w 



11/' V + 1) tia , - A/2 

W — - 1 — - — -W + re ' 



l -^V + l -H 0+ k 
r z 2 



(47) 
(48) 
(49) 



where a prime in these transformed equations indicates differentiation with respect to f. Note that there is no need to consider 
1391 1 because of the approximation X = 0. Besides, Eq. (I4H can be rewritten accordingly and manifest scale invariance. 

Similarly, we define the transformed Zerilli wave function Z — Z/M, which satisfies a scaled wave equation outside the 

star: 



d 2 



Vz(f.) 



Z(f») = 0, 



with the scaled potential Vz(f„) given by: 
1 - 2/f 



Vz(f.) 



[2n 2 (n + l)f 3 + 6n 2 f 2 + 18nr + 18] 



(50) 



(51) 



f 3 (nr + 3) 2 

It is clear that the scaled polar oscillation equations are explicitly mass-independent. Hence, the scaled polar QNM fre- 
quencies co of TVIIM could also display a similar universality as those of the axial oscillations, save cases in which perturbation 
in pressure becomes non-negligible, e.g. the pressure p-mode. Since the mass profiles of other realistic stars are close to that of 
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TVIIM, their QNMs should follow the same trend. In Fig. [5] and Fig. [3] the solid circles show the scaled frequencies of polar 
w-mode and /-mode QNMs of TVIIM, which clearly illustrate the scaling behaviour of these modes. Besides, they are also 
well approximated by the solid (dotted) line obtained from the best quadratic fit to QNMs of realistic neutron stars (TVIIM). 
The similarity between these lines in turn establishes independent corroboration of our theory. 

We have shown here that the scaled polar QNM frequencies Co of realistic neutron stars are approximately given by 
a universal function of the compactness. Our finding is particularly intere sting for the f-mode oscillatio n because Eq. JjJ 
is different from the empirical formulas and J1J discovered previously llAndersson fc Kokkotasl Il998l) . As mentioned in 
Sect. I, Eq. in general leads to explicit mass-dependence of Mlo t and is inconsistent with the scalin g behaviour of the polar 
oscilla tion equation. Moreover, comparing Fig. 3 in the present paper with Fig. 1 in the paper bv lAndersson fc Kokkotasl 
dl99S^ . we conclude that Eq. J5J can indeed provide a better fit to the scaled oscillation frequency Muj r . On the other hand, 
despite that Eq. @ is apparently different from @, it implies that Mlo\ = M/t is still a function of the compactness, namely: 

"«» (?) 4 Kf)+4 m 

Hence, it is again in perfect agreement with the scaling behaviour of the polar oscillation equation. It i s worthy of remark that 
the le ading dependence in (1521 . Mloi oc (M/R) i , can be argued from the quadrupole radiation formula jAndersson fc Kokkotasl 
1998) and is expected to be a good approximation as long as the compactnesss is small. In other words, Eq. 15211 is a Taylor 
expansion of Mun about the point M/R = 0, whereas Eq. @ expands Muj about a non-zero value of M/R (say, M/R « 0.2 
in the current situation). Therefore, the disparity between (j^J and 1521 1 is reconciled. As the compactness of realistic neutron 
stars usually lies in the range of 0.1 — 0.3, we expect that these two formulas can complement each other in the prediction of 
QNM frequencies of /-mode oscillations. 



6 APPLICATION 

As TVIIM successfully captures the universality in QNMs of w- and /-modes, it can be used to set up a benchmark for realistic 
stars. Bases on this, a viable scheme to infer the mass and the radius of a neutron star from its w- and /-modes QNMs is 
proposed here. In Fig.|Hl we plot 10 x Reojf/Reaj a (short dashed- line); 10 x ReoJf/ReoJp (long dashed-line) ; RetJa/Im^a (solid- 
line); and Reojp/Imojp (dot-dashed-line) against the compactness, where u a , io p and ut are QNM frequencies of the leading 
axial and polar w- modes, and /-mode of TVIIM, respectively. The curves are all monotonically increasing function of M/R 
in the regime 0.1 < M/R < 0.3 where physical neutron stars are stable. Depending on experimental data available, we can 
use any one of these ratios to approximately infer the compactness of an unknown neutron star. Once the compactness is 
known, we can readily read the approximate mass of the star from Figs. 1 — 5 if the real (or imaginary) part of relevant QNM 
frequency is obtained in gravitational wave detection. Since both the real and imaginary parts of the scaled frequencies of 
w- and /-modes are well approximated by quadratic functions of M/R, implementation of the above-mentioned procedure 
is straightforward. Still, feasibility of this scheme hinges on availability of accurate experimental data, which is beyond the 
scope of the present paper (see, e.g. iKokkotas et alJboOll for a discussion on this issue) . 

Here we illust rate our method with the following example. Consider a neutron star constructed from APR2 model 
l|Akmal et alJll998h . whose mass and radius are 1.596M0 and 11.77 km, respectively. The first polar w-mode QNM of this 
star has a angular frequency 64.57 kHz and a damping time of 2.666 x 10~ J s. Therefore, the ratio Reu> p /Imu> p is equal to 
1.722. If such wave signals are detected by us, we can use TVIIM as a benchmark for this star and infer its approximate mass 
and radius from the wave data. From the ratio Rewp/Imojp = 1.722 and the dot-dashed-line in Fig. [HI the compactness of the 
corresponding TVIIM star is found to be 0.2025 (in geometrized units). Hence, the scaled frequency can be obtained from 
Fig. |21 By comparing the scaled frequency with unsealed one, the mass of the star can be readily found, which is 1.581Mq in 
this case. Similarly, the radius of the star can be found from the compactness given above and is equal to 11.51 km. Both the 
values of the radius and the mass inferred from TVIIM star are close to those of APR2 star. From this example one can see 
that TVIIM can in fact serve as a template for other realistic stars. In this regard, the role of TVIIM in gravitational wave 
analysis is reminiscent of that of the hydrogen atom in atomic spectroscopy. 

Before we can put this scheme into practice, we have to further examine whether the universality demonstrated in 
Figs. 1 — 5 is truly generic. As discussed above, the universality is closely tied with the validity of TVIIM, which is developed 
under the assumption that the density function p(r) is expandabl e in series of r. Despite that the assumption is appealing, 
there is one well-known exception to it, namel y the strange star JWittenlll98l . One salient feature of the EOS of strange 
quark matter, say the MIT bag model (see, e.g. lGlendenninJll997l . chap. 8), is that the density is non-zero at zero pressure. 
Therefore, the density of a bare quark star (i.e. one without a normal nuclear matter shell) does not vanish at the surface 
of the star. Obviously, TVIIM fails to mimic such stars as there is a discontinuity in the density on the stellar surface. As 
shown in Fig. |2l the QNMs of bare strange quark stars deviate marked from the university demonstrated by other stars 
that are described by continuous density profiles. To pinpoint the cause of the deviation, the QNMs of strange quark stars 
with a normal nuclear matter shell are also shown in Fig. [0] It is obvious the QNMs of such hypothetically "dressed" quark 
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star, which have no discontinuities in their density profiles, also display the universality discussed in the present paper. We 
conclude that the scheme proposed above works properly as long as the star in consideration is characterized by a continuous 
density profile. Furthermore, by virtue of the TOV equations, the pressure is always a continuous function of r and the density 
profile is continuous if relevant EOS P(p) is a continuous function of p. As a result, deviation from the universal curve might 
indicate a discontinuous (or rapidly varying) EOS. This effect can indeed be observed from the QNM frequencies of GM24 
shown in Fig. [I] and Neutron stars constructed from GM24 will become unstable if M/R > 0.21 owing to emergence of 
phase transition that leads to discontinuous (or rapidly varying) EOS. Correspondingly, QNM frequencies of GM24 stars no 
longer follow the those of TVIIM as its compactness approaches that limit. 



7 CONCLUSION 

In this paper we have unveiled some universal scaling behaviours prevailing in QNM frequencies of w- and /-modes. In 
addition, we have shown that these behaviours are attributable to the generic form of the mass distribution function. Based 
on our discovery, a scheme is proposed here to infer the mass and the radius of an unknown star from its gravitational wave 
spectrum. Despite that there might be some exceptional cases, e.g. bare strange stars, the scheme does provide reasonable 
estimate of relevant values in normal situati on. As it has be en shown that the mass and the radius of a neutron star can 
be used to infer the EOS of stellar matter l|Lindblomlll99^ . our discussion will also lead to determination of EOS from 
gravitational wave signals. 

The TVIIM considered here is a two-parameter model that is characterized by the mass and the radius of a star. While 
reproducing the main features of realistic stars, it cannot capture the fine details in individual cases. Should it be possible to 
gather more than one QNM frequency of a pulsating star, we could, in principle, construct more refined template for it. We 
are currently working along this direction and relevant result will be reported in due course. 

Our analysis sketched above is based on assumption that QNMs of gravitational waves are clearly discernible. However, it 
is well known that gravitational wave sig nals are usually overwhelmed by noises of various kinds, including the rmal, quantum 
and gravity-gradient noises (see, e.g. lLiu fc Thorne!l200d iThorne fc Winsteinlll999llSantamore fc Levinlbool. and references 
therein). As yet there are no viable methods to completely eliminate these noises such that clear and reliable gravitational 
wave signals can be recorded. However, with the continuous efforts of researchers in this field, we will surely get rid of these 
noisy gremlins at a later date. Thus, gravitational wave astronomy will open up a new channel to observe distant compact 
stars and the scheme outlined here will be applicable in the near future. 
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Figure 1. (a) The real and (b) the imaginary parts of the scaled frequency Mtu of the least-damped axial ui-mode are plotted against 
the compactness M/R for neutron stars described by different EOS, including APR1, APR2, AU, GM24, Models A and C, UT, UU, and 
also TVIIM. The solid line represents the best quadratic fit in M/R to the scaled frequencies of the eight realistic stars. Likewise, the 
dottcd-line is the best quadratic fit to those of TVIIM. 
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Figure 2. (a) The real and (b) the imaginary parts of the scaled frequency Mu> of the least-damped polar ui-mode are plotted against 
the compactness M/R for neutron stars described by different EOS (see the caption of Fig. 1). The solid and dotted-lines represent the 
best quadratic fit in M/R to the scaled frequencies of the eight realistic stars and TVIIM, respectively. 
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Figure 3. (a) The real and (b) the imaginary parts of the scaled frequency Mu> of the fluid /-mode arc plotted against the compactness 
M/R for neutron stars described by different EOS (see the caption of Fig. 1). The solid and dottcd-lincs represent the best quadratic fit 
in M/R to the scaled frequencies of the eight realistic stars and TVIIM, respectively. 
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Figure 6. The mass distribution functions of neutron stars characterized by eight realistic different EOS and with a common compactness 
M/R = 0.2 are shown. The mass distribution function of TVIIM, the solid line, is shown as a benchmark for realistic stars. 
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Figure 7. The complex eigenfrequencies of the least damped polar M>-modc for stars constructed from seven realistic EOS (unfilled 
symbols) and TVIIM (solid-circle), which have a common compactness M/R = 0.2, are shown. The grey (or dark-grey for the case of 
TVIIM) symbols indicate the corresponding results obtained under the approximation X = 0. 




Figure 8. (a) lOx Rewf /Reoj a (short dashcd-line); (b) lOxRcajf/Retjp (long dashcd-line); (c) Rcu; a /ima; a (solid-line); and (d) ReUp/Imajp 
(dot-dashed-line) are plotted against the compactness M/R. 




Figure 9. (a) The real and (b) the imaginary parts of the scaled frequency Md of the least-damped polar ui-mode are plotted against 
the compactness M/R for neutron stars described by different EOS (see caption of Fig. 1). The solid and dottcd-lincs respectively show 
the results of "dressed" and bare strange quark stars. 



